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ABSTRACT 

Context. The TeV emitting X-ray binary LS I +61° 303 has two radio periodicities that correspond to a large periodic outburst with the 
same period as the orbit, 26.5 days (phase and a second periodicity of 1667 days (phase 0), which modulates the orbital phase and 
amplitude of the large outburst. Analyses of the radio spectral index revealed in LS I +61° 303 the presence of the critical transition 
typical for microquasars from optically thick emission (related to a steady jet) to an optically thin outburst (related to a transient jet), 
and found that it occurs at O „-,-,, which is modulated by 0: O c „, = /(0). 

Aims. We examine the possible implications of averaging high energy data over large and O intervals in the light of puzzling 
published INTEGRAL results, which differ for different averaging of the data. 

Methods. In microquasars, a simultaneous transition between two X-ray states occurs at the switch from optically thick radio emission 
to an optically thin radio outburst, from the low/hard to the steep power-law state. Assuming that the same transition occurs in 
LS I +61°303 at O cn „ we can show qualitatively the effect of averaging high energy data on ©, by analysing the effects of averaging 
radio spectral index data across the same © interval. We then model the two X-ray states, low/hard and steep power-law state, and 
show quantitatively how their mixing can affect the results. 

Results. When folded over too large a interval, spectral data from INTEGRAL can yield a false picture of the emission behaviour 
of the source along the orbit because it may be mixing two spectral states. Furthermore, averaging the data along the orbit may result 
in a dominant low/hard spectral state, which, for insufficiently extended sampling, might appear without a cut-off. 
Conclusions. The INTEGRAL results can be interpreted as two X-ray states that alternate with each other along the orbit. The 
implications of this analysis for HE/VHE data from LS I +61°303 and a possible connection between HE/VHE emission and the 
steep power-law state are discussed. 
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V) ■ 1. Introduction iMassi & Kaufman Bernadol d2009h analyzed the dependence 

of the radio spectral index a with flux density S k v", on both 

™ , ,. . ... v ,. T „ T £ , periods © and The analysis uses 6.7 years (1994-2000) of 

v N , The clear periodicity of the X-ray binary LS I +61 303, formed Jf TA „ . .™ T , „ , J T „ T , 10 -,;U . c 

— <• * . , .,/ ,/ „„,, , . _ NASA/NRAO GBI data from LS I +61 303 at two frequencies, 

by a compact object and a rapidly rotating B0 V star, is amongst 2 _ , „, „, „, , • , . 

./ . i- u • J L, ,.f , . , ' ° , vi = 2.2 GHz and vo = 8.3 GHz. The analysis had two main 



its most peculiar characteristics. The orbital period of P = 26.496 
days (phase <S) with Q> pe riastron - 0.23-0.28, Aragona et al. 2009; 
Casares et al. 2005) modulates the radio flux resulting in a large 
outburst towards apastron, whereas a superorbital period of 1667 



findings. The first was that in LS I +61°303 the optically thin 
outburst occurs after an interval of optically thick emission. In 
microquasars, the optically thin outburst is due to a transient 



i ^.uui^vn^^^u, "."TT""''"."."^ '"''^"^ 1 "" i et consisting of shocks propagating in a pre-existing optically 

r <j , days (phase 0) modulates both the amplitude and the orbital J tU . , .... ° . , . f„ , . , OAA ^irrr — nnS, , , 

-h J vr „ , , ,. ^ — I — nn7>r>nL 17^ 1 thick emitting steady let Fender et al. 2004; Massi 201 la). The 

03 ' phase of the large radio outburst (Gregory et al. 1999; Gregory , °. , . i. J ,. — f — f . 

- - ■ l onn-. k rr 5 1 ■ . c .u- — T — — 7 1 • — second remarka ble find ing resulting from the radio spectral m- 

2002). Two models exist for this system: the two-peak micro- , . „rrr — ■ c CV r ^ J^r>/-vrA ^ ^ • 

, , tt-z — ^ „ n , — i ir.r.k rr; zm — .1 l^,u*A dex analysis of Massi & Kaufman Bernado (2009) was that in 

quasar model (Marti & Paredes 1995; Bosch-Ramon et al. 2006; . „ T ,; ,„, , . . , , — — rr - ■ — , .77, , . w 

— — nr-irvrv4, — 7T — , — , , J,, , . i, r i LS I +61 303 this switch from one kind of a let (steady let) to 

Romero et al. 2007) and the pulsar model Maraschi & Treves , . . ... . , , . , 

inoi ilT-> 1. iMrvrv^ t-i c »u £ ! — , , the other one (transient let) occurs twice along the orbit at the or- 

1981; Dubus 2006). The name of the first model was coined be- .. . . j. , , J :, . . , , „ 

bital phases predicted by the two-peak microquasar model. One 



cause of the highly eccentric orbit in LS I +61 303, for which , • ,, , • 

J . , , . . , sees a spectral transition, optically thick/ optically thin emission, 

two ejections are expected along the orbit: one around penas- . r , . ' t, #>oT r.n jj 

J , , r , ,° . . _ „ r, . ■ . nrst around penastron at <P fn( 1 = 0.33 + 0.13 and then again 

tron and a second one shitted about AO ~ 0.3 towards apastron. , \ _ n n n ^ • . 1 

T1 ,. . .„ < < . ,j- V j towards apastron at <D C rit2 = 0.7 + 0.13. Owing to the qual- 

In the radio regime, the nrst model interprets the often observed - ' 

one-sided radio morphology and as well the observed switch to 

a double-sided morphology, as that of a precessing microquasar 

with a small and variable angle between the jet and the line of 

sight (Massi et al. 2004; Massi & Zimmermann 2010). The sec- 1 , , ■ 

, ... — ^7 ., , . . ., established relationship tor microquasars between radio and X 

ond model interprets the one-sided structures as a cometary tail ti r=, — " ^ TTZTUTii. 

of the pulsar (DhawaHFdD[2006). ray Stat6S ^nder et alj |200j). 



itative and quantitative agreement between radio data and the 
two peak microquasar model, the microquasar scenario is as- 
sumed here as a working hypothesis to explain the puzzling re- 
sults obtained from INTEGRAL observations. We use the well- 
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The two radio states, the steady jet (a >0) and the tran- 
sient jet (a <0), should be directly linked to two spectrally dis- 
tinct X-ray states in the unified model of X-ray states with ra- 
dio jets: the low/ha rd X-ray state and the steep power-law state 
dFender et al.l l2004). The X-ray states are clearly distinguishable 
by their energy spectra ([Grove et alj|1998t iFender et al. l l2004t 
iMcClintock & Remillardll2006l) . A steady jet (a > 0) occurs in 
the low/hard X-ray state and its X-ray spectrum is characterized 
by a power-law with photon index T * 1.5-1.8 and a cut-off 
at high energies. The high energy cut-off is different for differ- 
ent sources, with s ome showing a folding energy Ej a; 115 keV 
dGrove et al.lll998l) and others, as GX 339-4, a f olding energy 
Ef of only 66 keV (ICab allero-Ga rcia et al.ll2009T) . In particular 
for LS I +61°303, the very low X-ray luminosity, the observed 
optical modulation correlated with the radio outburst (see Massi 
& Kaufman Bernado 2 009), and the soft ening associated to a 
decrease in X-ray flux dSidoli et al.1 12006). which is consistent 
with accreting systems in quiescence (Corbel et al. 2006, 2008), 
all point towards the source exhibiting a very low hard, nearly 
quiescent, X-ray state. Finally, the transient jet (a < 0) is asso- 
ciated with the steep power-law X-ray state with a black-body 
component (with kT~l keV) and a power-law with a photon in- 
dex of F > 2.4 and no cut-off, thereby intrinsically allowing also 
for high energy emission. 

Data from LS I +61 "303 taken with the International 
Gamma-ray Astrophysics Laboratory (INTEGRAL/ISGRI 20- 
500 keV, LS I +61°303 is detectable in the 20-150 keV range), 
however, yielded surprising results, finding e.g. along most of 
the orbit a photon index only compatible with a low hard state, 
but without a cut-off. The photon index of the hard X-ray spec- 
trum, furthermore, depended on how much data was averaged 
toget her and whether it was also averaged over the orbital pe- 
riod ( Chernvakova et al.ll2006HZhang et al.ll2010l) . Here, we aim 
to explain these results by using the relationship between radio 
and X-ray states that allows us to predict the kind of emission 
expected at soft gamma-rays based on the radio emission, or 
viceversa. In particular, in LS I +61°303, the and i> period- 
icities allow us to test this relationship by using GBI data, even 
if no simultaneous observations are available (as is the case for 
INTEGRAL). From the radio spectral transitions seen in the GBI 
data one can infer which spectral transition is expected in X-rays 
and when. The combination of radio spectral index data from dif- 
ferent intervals can then help us to understand what happens 
when INTEGRAL data is averaged over large parts of 0. 

2. INTEGRAL observations versus and O 

In Fig. Q] INTEGRAL observations published to date are dis- 
played in terms of the intervals in which they were car- 
ried out (Zha ng et al.ll2010h and denoted as II, 12, and 13 (see 
also Fig. 1 in Zimmermann et al. 2011). Any spectral analysis 
should consider 0, because within the energy range covered by 
INTEGRAL data, two states, the low hard and the steep power- 
law state, can be detected, when they occur. 

Radio spectral index curves showing $> C rit,2, the transition 
for the large radio outburst towards apastron, for = 0.0-0.1 
and 0.7-0.8 are shown in Fig. [2] The orbital occurrence at phase 
$>crit,2 of the transition from a > to a < lies for = 0.0-0.1 
around w 0.73-0.77, and for = 0.7-0.8 it occurs around 
<b crit x 0.63-0.67. Owing to the connection between radio and X- 
ray states, it now becomes evident from Fig.|2]that, if the energy 
range of an instrument covers both the low hard and the steep 
power-law state, folding over the orbital period data of differ- 
ent 0s might result in the mixture of different states. The orbital 




1667 d phase (0) 



Fig. 1. Published high energy observations of LS I +61°303 by 
INTEGRAL (11,12 and 13) shown in the context of the 1667 day 
radio period (0). The underlying light curve gives 6.7 yr radio 
Green Bank Interferometer data at 8.3 GHz folded with 0. 



phase of the transition from optically thick emission to an opti- 
cally thin outburst does not remain the same over the 1667 d pe- 
riod and the unified X-ray state model with radio jets assumes a 
physical con nection between the radio spectral transition and the 
X-ray states dFender et al.ll2004T) . The folded data no longer pro- 
vides any accurate information about the orbital occurrence of 
the spectral transition. As we show here, the information about 
one state could even be lost completely, if e.g. the sampling for 
one interval is finer, which leads to the suppression of infor- 
mation from the other interval, and/or data from several different 
©-intervals, and not only the two intervals shown in Fig. [2] are 
mixed up. 

Folding all data or parts of II, 12, and 13 together in terms 
of orbital phase increases the sampling and the significance of 
the signal. In doing so, several authors established that the emis- 
sion in the i nterval 10-100 keV is clearly modulated by the or- 
bital phase d Chernvakova et al. 2 0061 iHermsen & Kuiperl l2007: 
IZhang et all I2010Q . Moreover. IChernvakova et all ([2006), who 
used data from almost all of 13 (covering = 0.57 - 0.99) 
and a few data points from the beginning of II (covering = 
0.0 - 0.06), have found that whereas along most of the orbit the 
photon index T was ss 1.4-1.7 (i.e. low hard state), it changed to 
F = 3.6^ J' j (i.e. compatible with a steep power-law state) in the 
orbital phase interval <I> = 0.6-0.8. This result is consistent with 
the ra dio spectral index analysis of Massi & Ka ufman Bernaddl 
(2009), giving for the same ©-intervals of 13 the transition to 
an optically thin outburst (a transient jet) after O = 0.6+0.1. 
In Fig. [2] one sees that for = 0.7 - 0.8 (which is part of 
13) the optically thin emission dominates in the interval <t> 
0.6 0.8 Following the radio spectral index though, folding over 
too much of a cycle could imply that, although the resulting 
light curves can establish the overall periodicity of the source at 
these energies, it m ixes different spectral states. The analysis of 
IZhang etall d2010l) . who used II, 1 2, and 13 together (f o llowin g 
the same data reduction method as IChernvakova et all ((2006)), 
thereby covering both the minimum and the maximum ampli- 
tude interval of 0, give a hard photon index of F ss 1.4-1.9 along 
the whole orbit (average: 1.7), no longer finding evidence of a 
very soft spectrum (F > 2.4) in the interval O = 0.6-0.8. 

How can these results be interpreted in terms of the radio 
spectral index analysis from above? By examining Fig. |2l we 
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Orbital phase 



Fig. 2. Radio spectral index for two differe nt intervals. The 
vertica l lines mark the <E> interval where IChernvakova et all 
(2006) determined a photon index compatible with a steep 
power-law state (see Sect. 2). 



can see, how by folding together data of two intervals, opti- 
cally thin and thick emission can be mixed, because the orbital 
occurrences of <b C rit,2, and therefore of the transition between 
optically thin/thick emission, differ for the intervals. The opti- 
cally thin emission then no longer dominates the orbital interval 
<t> = 0.6-0.8. When averaging over the intervals II, 12, and 13, 
optically thick and thin emission get mixed together in the <D 
interval of interest. Therefore, the measured photon index must 
not be directly interpreted as a continuous low/hard state of the 
source, because it represents the combination of different states. 
Moreover, we note that this mixing also occurs, when only one 
interval is selected and the data is averaged over the whole 
orbit. As one sees in Fig. [2] optically thick and thin emission 
occur alternately during the orbit. When averagin g the spectrum 
across the whole orbit, Chernvak ova et al.l d2006l) find again an 
average photon index of F ~ 1.6 but no cut-off at high energies. 
An example of the effect of this mixing of X-ray states is shown 
in Fig. |3j here we show one spectrum for a low/hard state (open 
squares) and one spectrum for a steep power-law state (open tri- 
angles) and the spectrum for the low/hard state is 



/(£) = 



E r exp -(£ 



E T E <E C 
' E^jEfoid E > E c 



(1) 



where E c is the cut-off ener gy and E f u the folding energy, 
which is in iGrove et"aT1 d 1998b assumed to be approximately 
twice E c . We use as an example the source GRS 1716-249 
with Efoid = 115 + 8 keV, because this source, analysed by 
IGrove et"a"fl ([1998) in both spectral states, results in a steep 
power-law state of significantly lower luminosity. This might 
also hold for LS I +61° 303, as suggested by the results of 
IChern vakova et al. (2006) (i.e. dominance of the low hard state 
when mixing the data). In Fig. [3] we superimpose an average of 
the two curves (filled circles). The four vertical lines indicate the 
fo ur points of the averaged (o ver both O and 0) spectra derived 
bv IChernvakova et al.l (|2006). The fit to these points resulted in 
a photon index F consistent with a low hard state, but no cut-off 
at high energies was found. This example shows how an average 
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Low hard state with r=1 .53 

+ cut-off (E-fold~1 15 keV) 

Steep power law with r=2.42 

+ Blackbody component (kT=1 keV) 

Average (r=1 .57 for 20-80 keV and r=1 .53 

for 0.1-20 keV) 




0,01 
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Fig. 3. Model spectra for the low/hard state (power-law plus cut- 
off, open squares), the steep power-law state (open triangles) to- 
gether with a black-body disk component with kT~ 1 keV and 
the average of the two states (filled circles). 



could yield a curve, where the photon index of the power-law 
below the cut-off energy (E c « 115/2 keV) is consistent with a 
low/hard state although the steep power-law is part of the data. 
In addition, we see in this example for GRS 1716-249 that onl y 
the last point, among the four ones bv lChernvakova et alJ (l2006). 
samples the average curve after the cut-off and lies, with around 
84 keV, still quite close to the cut-off energy (E c « 115/2 keV). 
It would be difficult under this condition to detect the cut-off. For 
GRS 1716-249, Grove et al. indeed established that there is ex- 
ponential decay with Ef id = 115 + 8keV for very well-sampled 
data up to 300 keV. However, if the cut-off in LS I +6 1 °303 were 
to occur well before the assumed one, the contribution of the 
steep power-law would strongly bias the analysis of the averaged 
data. 



3. Conclusions and discussion 

We have investigated the implications of the radio spectral in- 
dex transition, O cn -, = /(0), in LS I +6T303 for the analy- 
sis of INTEGRAL data. In microquasars, this transition is di- 
rectly linked to the transition between two spectrally distinct 
high energy st ates that cover t he energy range from X-rays to 
gamma-rays (IGrove et alj 19981) . but which in this paper are sim- 
ply referred to as X-ray states: the low hard X-ray state and the 
steep power-law X-ray state. For the energy range covered by 
INTEGRAL data, both states can be detected, when they occur. 
From our analysis, we conclude that: 

- Folding data over too large intervals could mix up data 
from different states so that the resulting spectra can no 
longer be unambiguously interpreted and mimic a persistent 
low hard state. 

- For data from the energy range accessible to INTEGRAL, 
even if the data were not folded on 0, averaging across <D 
will destroy any information about the different states. 

- We have pointed out that the lack of a cut-off in the power- 
law for data averaged all over the orbit, means that if 
LS I +61°303 has a cut-off well below -115/2 keV, the con- 
tribution of the steep power-law will strongly bias the analy- 
sis of the spectrum of the averaged data. However, we have 
shown that a higher cut-off requires highly sensitive sam- 
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pling, extending well beyond the used range of up to 100 
keV. 

The INTEGRAL results are then explainable with two X-ray 
states alternating along the orbit in parallel to, as expected in 
the microquasar scenario, the two radio states. What does this 
imply for the analysis of the HE/VHE data? For VHE, when 
LS I +6T303 is detected by MAGIC or VERITAS, its spectrum 
is always closely fitted with a photon index of >2 .4 ( Albert "et alj 
2009; lAcciari et al.l 120111: IMAGIC Co llaborati on et alJ 1201 lh . 



which is the same as expected for the steep power- law. In 
this context, the observations of Anderhub TaTT(l2009T) are of 
special interest. They have observed simultaneously X-ray and 
VHE emission from LS I + 61°303 and found a correlation. 
M assi & Zimmermannl d2010l) noticed that for these observa- 
tions the two fluxes result in Fvhe k w ^ tn f = 0.99, 
which is in agreement with t he co rrelation observed in blazars 
dKatarzvnski & Walczewskal l2Q10t> . Here the X-ray emission 
is due to synchrotron, and VH E is synchrotron self-Compton 
(SSC) dKatarzvnski et alj|2005l) . For LS I +61°303, several au- 
thors have explained the X-ray excess around apastron with 
synchrotron emission and the VHE wi th either external inverse 
Compton (EIC) or SSC emission dGupta & Boettcherl [2006; 
IZabalza et aT]|201 llhlZabalza et all d201 lh demonstrated that the 
X-ray/VHE correlation observed by Anderhub et al. (2009) is 
compatible with a one-zone leptonic particle population produc- 
ing the emission. Their models use the observed X-ray photon 
index (1.55-1.67) and can reproduce the observed lightcurves 
quite well. Interestingly, this X-ray photon index is comparable 
to the one expected for the low hard state and could be mistaken 
as such. The soft X-ray synchrotron emission during the steep 
power-law state can then explain two findings. First, the VHE 
emission by invoking EIC or SSC and, second, the absence of 
an accretion disk black-body component in the soft X-rays (0.1- 
10 keV) expected during the steep power-law state, if the syn- 
chrotron X-ray component outshines the disk. Owing to the lack 
of simultaneous hard X-ray/TeV observations of LS I +61°303, 
it remains open, if the jet component extends into the hard X-ray 
regime. 

For the high energy emission in the GeV r ange, detected 
e.g. w ith Fermi-LAH, the picture is more diverse. IZabalza et alj 
d201 lh note that in their one-zone population model, a broken 
power-law electron energy distribution is needed to explain the 
GeV component as observed with Fermi-LPR. The steep power- 
law might not then be the only emission process, also because 
GeV emission is detected all along the orbit. Electrons from the 
steady jet can always upscatter stellar UV photons to GeV ener- 
gies (i.e. EIC see Bosch-Ramon et al. 2006). Nevertheless, more 
energetic particles from the transient jet could in addition pro- 
duce as well GeV emission via EIC and SSC. Intriguingly, the 
spectrum measured by FermZ-LAT shows, in addition to a power- 
law with a cut-off around 6 GeV, upper limits possibly compat- 
ible with the spectrum measured with MAGIC and VERITAS 
(see e.g. Fig. 2 in Hadasch 201 1). 

In any case, owing to the connection between the radio 
and the HE/VHE emission, the analysis performed here for 
INTEGRAL can and should in principle be extended to other 
high energy instruments, if the © coverage is good enough. In 
particular, the flux evolutions recently detected at HE/VHE im- 
ply that this data should b e analysed on as short a i nterva l 
as possible. As discussed in Massi & Kauf man Bernaddl (|2009), 
the two peak shape of the a versus <D curve varies with ©. 
It is most pronounced in ©=0.0-0.1 and 0.9-1.0, appreciable 
at 0=0.2-0.3 and 0.7-0.9, and absent at 0=0.3-0.6 (see Fig. 5 



in Massi & Kaufman Bernadd 2009). In addition, the distance 
of the two peaks varies with 0. These variations in the accre- 
tion curve are pr e dicted by the two-peak microquasar model of 
Marti & Paredes (1995) by incorporating variable wind veloc- 
ities of the Be star (see their Fig. 6). A reflection of this be- 
haviour at other wavelengths owing to the related emission pro- 
cesses is expected. As a matter of fact, there is an interesting 
increase in the overall flux level observed with Fermi-LKT af- 
ter March 2 009 (0 & 0.92 ) together with a broadening of the 
peak shape (Hadasch 2011). The flux variations observed with 
Fermi-LAI to date are therefore consistent with a varying with 
0. Similarly, strong variations are observed at very high ener- 
gies. The source went from being detected around apastron with 
VERITAS and MAGIC (VERITAS: 0=0.38, 0.4-0.47; MAGIC: 
0=0.17-0.25, 0.37-0.44, 0.59-0.6) to becoming quiescent dur- 
ing 2008-2010. No detection was reported by VERITAS for 
0=0.83-0.08, but weak detection was detected by MAGIC for 
0=0.05-0.1 1 around apastron. The source was again detected by 
VERITAS in October 2 010 corresponding to 0=0.26, but this 
time around periastron (Albert et al.l 120091 ; lAcciari et alj|201 it 
IMAGIC C ollaborati on et alj|201 ll) . Nonetheless, at the moment 
the insufficient coverage of these instruments does not allow 
a closer comparison. Nevertheless, the observations clearly indi- 
cate that for the analysis of these data, should be considered 
to avoid any mixing of signal from the outbursts. 
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